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Abstract 

Ultrafast time-resolved x-ray scattering is an emerging approach to image the dynamical evolu- 
tion of the electronic charge distribution during complex chemical and biological processes in real- 
space and real-time. Recently, the differences between semiclassical and quantum-electrodynamical 
(QED) theory of light-matter interaction for scattering of ultrashort x-ray pulses from the electronic 
wavepacket were formally demonstrated and visually illustrated by scattering patterns calculated 
for an electronic wavepacket in atomic hydrogen [Proc. Natl. Acad. Sci. U.S.A., 109, 11636 
(2012)]. In this work, we present a detailed analysis of time-resolved x-ray scattering from a sam- 
ple containing a mixture of non-stationary and stationary electrons within both the theories. In 
a many-electron system, the role of scattering interference between a non-stationary and several 
stationary electrons to the total scattering signal is investigated. In general, QED and semiclas- 
sical theory provide different results for the contribution from the scattering interference, which 
depends on the energy resolution of the detector and the x-ray pulse duration. The present findings 
are demonstrated by means of a numerical example of x-ray TRI for an electronic wavepacket in 
helium. It is shown that the time-dependent scattering interference vanishes within semiclassical 
theory and the corresponding patterns are dominated by the scattering contribution from the time- 
independent interference, whereas the time-dependent scattering interference contribution does not 
vanish in the QED theory and the patterns are dominated by the scattering contribution from the 
non-stationary electron scattering. 
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I. INTRODUCTION 



To fully understand the functionality and dynamic behavior of molecules, solids and 
complex biological systems, it is important to image the motion of electrons in real-time 
and in real-space. The motion of atoms within molecules and solids that is associated 
with chemical transformations occurs on the femtosecond (1 fs = 1CT 15 s) timescale. The 
timescale of electronic motion, responsible for electron-hole dynamics and electron transfer 
processes in molecules can be even faster, on the order of attoseconds (1 as = 10~ 18 s) lj- 
5[. The ultimate goal of the emerging field of time- resolved imaging (TRI) is to visualize 
electronic motion on an ultrafast timescale as electrons move in atoms, complex molecules or 
solids, as occurring for instance in photoinduced exciton dynamics, during bond formation 
and breakage, conformational changes and charge migration 6r|9|. With the tremendous 
advancement in technology for producing ultraintense and ultrashort x-ray pulses from novel 
light sources, it seems possible to obtain information about ultrafast dynamics of electrons. 
In extension of the concept of "molecular movies", which track the motion of atoms on fs 



timescale 
(FEL) 



13] . ultrashort, tunable, and high-energy x-ray pulses 



15|, laser plasmas [16] and high-harmonic generation 
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rom free-electron lasers 



181 ] promise to provide 



"electronic movies" that take place on few fs to as timescale 19] . Recent breakthroughs make 



it possible to generate hard x-ray pulses of a few fs 



0,0, 



and pulse duration of 100 as can in 



principle be realized 2ll l22|. Utilizing the remarkable properties of x-rays from FEL, several 



new insi ghts have been gained about systems ranging from atoms 23j, |24j, molecules 25l. l26|. 
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29[, to matter in extreme conditions 30J, |31 |. 



clusters 127], complex biomolecules 

n 

Since the discovery of x-rays [32| , scattering of x-rays from matter has been used to unveil 



the structure of molecules, solids and biomolecules with atomic-scale spatial resolution 
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37J. Also, scattering of x-rays from atoms and molecules has been pro pos ed to gain insight 



38 



In order to 



about the excited electronic states of atomic and molecular systems 
image electronic motion in real-time and real-space with spatial and temporal resolutions of 
order 1 A and 1 fs, respectively, one can perform scattering of ultrashort x-ray pulses from 
the dynamically evolving electronic charge distribution. Pump-probe experiments are the 
most direct approach, where first a pump pulse induces the dynamics and then subsequently 
a probe pulse interrogates such induced dynamics. By varying the pump-probe time delay, 
one obtains a series of scattering patterns that serve to image the electronic motion with 
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atomic-scale spatio-temporal resolution. 

Recently, the theory of TRI from a non- stationary electronic system, using both the semi- 
classical and the quantum electrodynamical (QED) treatment of light-matter interaction, 



has been developed 



40] . In a semiclassical theory of light-matter interaction, matter is 



treated quantum mechanically and light is treated classically. In such a situation, the time- 
dependent Schrodinger equation is solved for the electrons together with Maxwell's equations 
for the light. By solving Maxwell's equations for given charge and current densities, the ex- 
pression for the differential scattering probability (DSP) is obtained, which is a key quantity 
in x-ray scattering. According to the semiclassical theory, x-ray TRI would be expected 
to provide access to the instantaneous electron density of the non- stationary electronic sys- 
tem. On the other hand, in a consistent quantum theory of light-matter interaction, where 
both matter and light are treated quantum mechanically, x-ray TRI encodes the information 
about spatio-temporal density-density correlation. Both the theories have been applied to 
an electronic wavepacket prepared as a coherent superposition of eigenstates of atomic hy- 
drogen and it has been shown that the scattering patterns obtained using both the theories 
differ drastically from each other. Moreover, it was shown that the patterns obtained within 
QED theory follow the motion of the wavepacket providing the correct periodicity of the 
motion, which cannot be captured by the semiclassical theory. In that case, the notion of 
the instantaneous electron density as the key quantity being probed in x-ray TRI for a suffi- 
ciently short pulse completely breaks down 40|]. However, x-ray TRI from a single isolated 
hydrogen atom is not a realistic scenario. In practice, a sample contains several electrons 
and when a tunable pump pulse with broad energy bandwidth interacts with an iV-electron 
system, one or few electrons participate in the formation of an electronic wavepacket and 
other electrons remain stationary. In such a situation, when an x-ray pulse scatters from 
a sample containing one or more non-stationary electrons and several stationary electrons, 
there is no way to know whether the scattering has taken place from the non-stationary elec- 
trons or from the stationary electrons and how the two scattering paths interfere with each 
other, i.e., interference between scattering from non-stationary and stationary electrons in 
the scattering process. Therefore, at this juncture it is important to analyze different types 
of contributions to the total scattering signal. The total signal can be decomposed into three 
main parts: first from stationary electrons, second from non-stationary electrons and third 
from the interference between non-stationary and stationary electrons. In the present work, 
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we will analyze how these different contributions in an iV-electron system contribute to the 
total scattering signal in both the theories (semiclassical and QED). 

This paper is structured as follows. Section II discusses the formalism and results for 
x-ray TRI in the case of many-electron systems, where only one electron forms an electronic 
wavepacket and other electrons serve as stationary reference scatterers in both the theories. 
Effects of different parameters such as energy resolution of the detector, pulse duration 
and spectral bandwidth of the x-ray pulse etc. for the electron-electron interference in 
the scattering process are discussed in detail. Section III presents a numerical example 
of x-ray TRI for an electronic wavepacket in helium, where one electron forms a coherent 
superposition of one-electron eigenstates and the other electron remains stationary and serves 
as a reference scatterer. In this particular situation for helium, the role of the scattering 
interference is investigated. Conclusions and future outlook are presented in Sec. IV. 



II. THEORY 



Our investigations are based on the theory for x-ray TRI of electronic wavepacket motion 



as developed in Ref. 



40|. Our equations are expressed in atomic units |4l|. Under the 



assumptions that the probe pulse is centered at the energy of the incident pulse with very 
small energy width and the coherence length of the pulse is large in comparison to the size 
of the object, the expression for the DSP within the semiclassical theory is related to the 
Fourier transform of the instantaneous electron density, p e (x, t), as follows 40 1 

2 

(1) 



dP 



dP P 



J d?x /9 e (x, t) e* Q ' x 



Here, 4& is the DSP for a free electron and Q is the photon momentum transfer. Here, 
we have assumed that the x-ray pulse duration is shorter than the dynamical timescale of 
the electronic wavepacket. According to Eq. (JT]), the measured scattering pattern provides 
access to the instantaneous electron density as a function of the pump-probe delay time t. 

Let us consider a scenario for time-resolved scattering of ultrashort x-rays in order to 
image the motion of a one-electron wavepacket in the presence of N stationary electrons. 
In such a situation, the N stationary electrons serve as reference scatterers in the total 



scattering signal. In this case, using the language of second quantization 



42|, the total 
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electronic wavepacket can be written as 



|vI/(t))=c; (4) |$o) = ^a a e-^ct|$ ), 



(2) 



with 



|$o> = j J> J |$o) 



(3) 



Here, cj, (c p ) creates (annihilates) an electron in spin orbital \ip p ) and e p is the orbital energy 
corresponding to \<p p ), i.e., H\<p p ) = e p \ip p ). H represents the electronic Hamiltonian at the 
mean-field level and |$o) is the unperturbed ground state of the iV-electron system with the 
electrons filled to the Fermi level. Here and in the following, indices p, q,r,s,... are used for 
general spin orbitals (occupied or unoccupied). Occupied orbitals in |$ ) are presented by 
indices i,j,k,l,..., whereas unoccupied (virtual) orbitals are symbolized by a,b,c,d, . . .. 
We rewrite the key quantity in Eq. ([1]) in terms of the density operator 

2 



d 3 x / rfV(^(t)|n(x')|*(t))(^(t)|n(x)|^(t))e lQ ' 



iQ-(x-x') 



with 
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x) = ^(xMx) = ^y?J(xVg(x; 



CpCg. 



(4) 



(5) 



pq 



Here, h(x) is the electron density operator, and the field operator tp^(x) [^(x)] creates 
(annihilates) an electron at position x. Using the expression for the wavepacket as introduced 
in Eq. (jSJ), Eq. (J4]) simplifies as follows 

2 

(6a) 



J d 3 xp e (x,t) e iQ x 



i 

+ E E <a b e^-^{q t £ ab + C u C* ab } 



i ab 



J2<u b £abe l{£a - £b)t 



ab 



with 



and 



I"! 



pq 



<Pq\ 



rf 3 x^(x)e- jQ -> 9 (x). 



(6b) 
(6c) 

(7) 
(8) 



6 



Here, the right-hand side of Eq. (I6al) provides the time-independent contributions due to 
scattering from the N stationary electrons. The second term as shown in Eq. f!6bj) is due 
to the scattering interference between the N stationary electrons and the non-stationary 
electron. The last time-dependent term in Eq. f )Bc|) is solely due to scattering from the 
non-stationary electron. 

On the other hand, in the full quantum theory of x-ray TRI, both matter and x-ray 
pulse are treated quantum mechanically and first-order time-dependent perturbation theory 
is employed for the interaction between matter and x rays. Here we assume that the probe 
pulse has a small bandwidth and a small angular spread so that the pixel assignment is well 
defined in the momentum space, the coherence length of the pulse is large in comparison 
to the size of the object, and the pulse duration should be sufficiently short to freeze the 
dynamics of the electronic wavepacket. Under these assumptions, the resulting expression 
for the DSP from a coherent, Gaussian x-ray pulse is j^O] 

dp d Pe r , t t 7 / \ r dr „, x _ i( )T 



poo . . poo J 

d£l dVl J w km 2tt 



x / (fx / dV ( * [t + - n (x ) e- inr n (x) 




e iQ.(x-x')_ (g) 



Here, C(r) = exp(— 21n ^ T2 ) i s a function of the pulse duration tj. tOk in and Wk s refer to 
the energy of the incident and scattered photon, respectively, while WajsC^O is a spectral 
window function centered at u^ in with a width AE. WaeC^O m °dels the range of energies 
of the scattered photons accepted by the detector. 

In the case of QED theory, an energy-resolved scattering process is considered for x-ray 
TRI. Therefore, any inelastic (Compton) scattering contributions due to excitations from the 
iV stationary electrons can be easily distinguished by utilizing the energy-resolving detector, 
if we assume that AE is small in comparison to the characteristic excitation energies of 
|$o)- Hence, excitations from the iV stationary electrons are not considered in the following. 
Similarly, on using the expression for the wavepacket as introduced in Eq. (j2J), the key 
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expression of Eq. fl9]) is simplified as follows 



d 6 x I d 3 x' (Wit 



T 



h fx') e~ iHT h fx) 



# it 



3 iQ.(x-x') 



06 



(10a) 
(10b) 
(10c) 
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Here, the first term, Eq. ( llOaj) . provides the time- independent contribution due to scatter- 
ing from the iV stationary electrons, which is identical to Eq. (l6aj) . The second term in 
Eq. ( 110b!) is due to the scattering interference between the N stationary electrons and the 
non-stationary electron, which seems different to the one shown in Eq. ( I6bl) . The last time- 
dependent term, Eq. fllOcj) . is again solely due to scattering from the non-stationary electron. 
The scattering contributions from the non-stationary electron, Eqs. fl6c|) and f llOcp . are not 
identical and provide completely different information about the electronic motion as shown 
in the case of a one-electron wavepacket in atomic hydrogen 40(. It is evident from Eqs. (J6j) 
and ( flOl) . that the first term in both the theories, Eqs. (l6al) and (llOaj) . provides identical 
scattering contributions to the total signal. In many systems of interest, the number of 
stationary electrons is large and therefore, the time-independent terms contribute a strong 
static background in the total signal. Due to the large number of the stationary electrons 
and one or few non-stationary electrons, the dominating time-dependent contributions in 
the total scattering signal are due to the scattering interference between stationary and 
non-stationary electrons (unless one is considering Q for which Yli^-a * s small). Therefore, 
it is crucial to analyze the scattering interference term, Eqs. fl6b[) and fll0b|) . in both the 
theories. 

On substituting Eqs. f!6bj) and fllObj) in Eqs. ([1]) and ()9]), respectively, the expression for 
the scattering interference contribution to the DSP, ^§fS can be written as 

dP mt dP e ^ ^ a:ab Se a -e b)t{clCab + ^jr^ (U) 



dtt 



ab 
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in the semiclassical theory, and in the QED formalism can be expressed as 

{2 _ 2 2 _ 2 ^ 

£J i £ a6 e-^( Wk -- Wk in +(5 V t) ) +A 4 £: fe e-^K-^- (£ ^ ) ) I. (12) 

It is important to analyze several aspects of both expressions, Eqs. fill I) and (j!2p . under 
different circumstances and properties of the probe pulse. 

• The Gaussian distributions in Eq. (1T2"]) are centered at positions u-^ in + ( Sa ~ Sb ) and 
<^ki„ — ( £a ~ £ '' ), and the width of the distributions is determined by tj. Therefore, 
if A.E is larger than the separation between the two distributions and the detector 
has a lower energy resolution than the energy spectral bandwidth of the x-ray pulse, 
then all the scattered photons are detected by the detector and Wae is constant in 
the energy range contributing to the integral in Eq. ( fi~2"l) . Hence, on performing the 
energy integral, Eq. ( TT2l) reduces to Eq. ( ITTj) . However, when the energy resolution 
of the detector is poor, it is very difficult to assign a unique pixel in Q-space to the 
scattered photon in the detector and it may no longer be possible to filter out Compton 
scattering from the N stationary electrons. 

• If the probe pulse is very short in comparison to the dynamical timescale of the elec- 
tronic motion, the r dependent exponent in Eq. ( 110b!) will reduce to unity. Therefore, 
if AE is smaller than the energy spectral width of the pulse and centered at the in- 
cident energy such that o>k s ~ w k in , the energy integral, Eq. ( TT2l) reduces to Eq. ( ITT1) 
and the scattering interference contributions from both the theories will be identical. 
However, an ultrashort pulse corresponds to a large energy spectral width due to the 
energy-time uncertainty relation. Therefore, the unavoidable energy bandwidth of the 
pulse causes an uncertainty in the momentum distribution of the incoming photon. In 
such a situation uniqueness of a pixel in Q-space is lost. 

• In general, however, the probe pulse is not very short in comparison to the dynamical 
timescale of the electronic motion. In such situations the r dependent exponent in 
Eq. fllObj) will not reduce to unity and both expressions for the interference will provide 



different contributions to the total signal. 



9 



In the following subsection, we present an example of one-electron wavepacket motion in 
helium. In this example, we apply both the approaches to compute time-resolved scattering 
patterns and analyze the contributions from the scattering interference to the total scattering 
signal. 



III. NUMERICAL RESULTS AND DISCUSSION 

A schematic scenario for probing an electronic wavepacket motion in helium is shown in 
Fig. [TJ The ground state configuration for both the electrons is Is 2 . The ionization potential 
of the first electron is 24.59 eV and 54.42 eV for the second electron 43] . A pump pulse with 



broad energy bandwidth excites one of the electrons from the ground state configuration and 
prepares a coherent superposition of the ls3d and ls4f configurations with the projection 
of orbital angular momentum being equal to zero. The energy difference between the Is 2 
and ls3d configurations is 23.07 eV and between the ls3d and ls4f configurations is 0.66 
eV (see Fig. [T]). Therefore, the dynamical timescale of the electronic wavepacket motion, 
which is inversely related to the energy spacing between the eigenstates participating in 
the wavepacket, is 6.25 fs. The spatial extension of the wavepacket is 14-17 A along the z 
axis and 7.5-9 A along the x and y axes. It is known that for orbital angular momentum 
quantum number equal to or larger than two, the quantum defect is almost zero. Therefore, 
the electron in the Is orbital sees no shielding of the nuclear charge and the electron in the 
superposition of 3d and 4/ orbitals sees complete shielding of the nuclear charge. Therefore, 
the wavefunction for the stationary electron in the Is orbital can be expressed in terms of the 
hydrogenic wavefunction with nuclear charge Z = 2. On the other hand, the wavefunction 
for the non-stationary electron in 3d and 4/ orbitals as well as other higher-lying orbitals 
can be written in terms of the corresponding hydrogenic wavefunction with Z = 1. This 
type of procedure for treating the two-electron problem has already been used in the past 
and was successfully applied to describe different types of physical processes |44j. 



In order to compute time-dependent scattering patterns of the electronic wavepacket (cf. 
Fig. [1]) as a function of the delay time in helium, we employ both the semiclassical and the 
QED approaches, i.e., Eqs. ([1]) and (J9]). Since the non-stationary and stationary electrons 
are energetically distinguishable (energy difference is around 23 eV) and an energy-resolved 
scattering process is considered with energy resolution at least equal to the unavoidable 



10 



T 




FIG. 1: A schematic scenario for probing an electronic wavepacket motion using ultrafast time- 
resolved x-ray scattering in helium atom. A pump pulse with broad energy bandwidth (indicated 
in orange) excites one of the electrons from the ground state and forms a coherent superposition 
of the ls3d and ls4f eigenstates with equal population. An isosurface of the electronic charge 
distribution of the wavepacket is shown (indicated in red), which undergoes periodic oscillation 
with oscillation period T = 6.25 fs. The dynamically evolving electronic charge distribution is 
probed by an ultrafast x-ray pulse (indicated in blue) . By varying the pump-probe time delay, one 
obtains a series of scattering patterns that serve to image the electronic motion with atomic-scale 
spatio-temporal resolution. 

spectral bandwidth of the probe pulse, any excitation from the stationary electron can be 
easily filtered out and therefore is not considered in the present case. The scattering operator 
e «Q x j g eX p anc l ec } j n terms of the spherical Bessel functions, ji(Qr), and spherical harmonics 

YM(f>) as 

e*Q* =4tt^ 3l (Qr) l£*(a,/9) (13) 

l,m 

where r = |x|. After introducing hydrogenic wavefunctions, the expression in Eq. (Q fac- 
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torizes into radial and angular parts. The angular part is given by 

h hh\( h h h ! _ (14) 

I 1 mi m 2 m 3 

whereas the radial part is numerically integrated. To calculate the patterns as a function 
of the delay times, we used a Gaussian pulse of duration 1 fs with 4 keV incoming photon 
energy, and assumed a Gaussian photon energy detection window of width AE = 0.5 eV 
for the detector. The patterns are calculated for Q ma , x = 2 A -1 corresponding to a 3.14 
A spatial resolution and to a detection angle of scattered photons of up to 60°. In order to 
compute the scattering contribution from the non-stationary electron to the total scattering 
signal, all transitions induced during the scattering process within the energy detection 



window are computed 40] . Therefore, transition amplitudes from the eigenstates involved 
in the electronic wavepacket to all the electronic states within the detection range of AE 
are computed, which includes all types of multipole transitions allowed by the conservation 
of angular momentum and amounts to 22000 transition amplitudes. Also, the scattering 
patterns using semiclassical theory are calculated by convolution of the square of the Fourier 
transform of p e ( x j t) with an x-ray pulse of duration 1 fs. 

Scattering patterns in the Q x - Q z plane (Q y = 0) as a function of the delay time at 
times 0, T/4, T/2, 3T/4, and T are depicted in Fig. [2j The time-dependent patterns 
shown in Fig. [2]^a) are computed within QED theory using Eq. ([9]), whereas in Fig. [2]^b) are 
computed within semiclassical theory using Eq. (CQ). It is evident from Fig. [2]^a) that the 
patterns undergo spatial oscillation along Q z in momentum space and reflect the motion of 
the wavepacket along z in real space (see Fig. [1]). One of the striking features of the patterns 
shown in Fig. |2]^a) is that when the charge distributions are symmetric, and corresponding 
patterns are asymmetric and vice versa, which can be understood as follows. The charge 
distributions are identical at delay times T/4 and 3T/4, as may be seen in Fig. [TJ while the 
electron clouds move in opposite directions at the two times. At time T/4, the flow of the 
electron cloud is downwards, whereas at time 3T/4 the flow is upwards. This is reflected 
by their corresponding patterns. Therefore, the patterns calculated within the QED theory 
capture the dynamics of the momentum distribution of the wavepacket. As a consequence, 
the apparent motions of the charge distributions and of the scattering patterns are shifted 
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FIG. 2: Scattering patterns in the Q x - Q z plane (Q y = 0) of helium (cf. Fig. 1). Scattering 
patterns obtained using (a) QED theory, i.e., Eq. ([9]) and (b) semiclassical theory, i.e., Eq. ([I]) at 
pump-probe delay times 0, T/4, T/2, 3T/4, and T, where the oscillation period of the electronic 
wavepacket is T = 6.25 fs. 

by 90°. On the other hand, the patterns shown in Fig. EJ^b) do not change significantly as a 
function of the delay time. 

In order to understand the scattering from the non-stationary electron, the scattering 
patterns corresponding to the non-stationary electron are shown in Fig. [3] in the Q x - Q z 
plane (Q y = 0) as a function of the delay time. The patterns shown in Fig. |3Ja) are obtained 
using Eq. f llOcp and in Fig. |3jb) are obtained using Eq. ( I6c|) . It is evident from Fig. [3ja) that 
the scattering patterns undergo oscillations as a function of the delay time. The scattering 
patterns shown in Fig. [3jb) are localized in the low Q region, which reflects the spatial 
extension of the electronic charge distribution of the wavepacket and also undergo changes 
as a function of the delay time, but do not display oscillations as the charge distribution 
oscillates. Hence, patterns obtained using semiclassical theory provide half of the actual 
period of the motion as the patterns start repeating themselves in half of the actual time of 
the motion. 
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We now investigate the contribution from the scattering interference between the sta- 
tionary and non-stationary electrons to the total scattering signal in both the theories. On 
re-writing Eqs. fl7|) and (|S]) in the case of helium, we find that 

C U u = J rf 3 Vi s (x)e jQ > ls (x) (15) 
is a purely real number, as is evident from Eqs. (|T3|) and (I14p . whereas 

C 3dAf = J rf 3 x^(x)e jQ > 4 /(x) (16) 
is a purely imaginary number, which also follows from Eqs. (|T3|) and (|T4|) . Therefore, 

and 

In a similar way we can write £3^ = £3^, an d £4/4/ = ^4/4/ > which are both real. On 
substituting the contributions of the £'s in Eq. fl6bl) . the contribution from the scattering 
interference within the semiclassical theory can be written as 

E <^-^{cl sls C ab + C lsls Cl b ), (19) 

a=3d,4/ 6=3(2,4/ 

which can be further decomposed into two parts: a time-dependent scattering interference 
contribution 

«3d«4/e i(e --^ )4 {£Li^3d4/ + Aa^aW) + «l/«3de J ^- e -)'{£* sls £ 4/3d + £ lsls £: /3d } 
= a^fe^-^Cuu^Mf - £ 3d4 /} + al/Q^e^^'A^-^^ - A/m} = 0, 

which is zero, and a time-independent scattering interference contribution 

2£i s i s [|a 3( i| 2 £3d3<i + I a 4 / 1 2 £4/4/]. (21) 



Similarly on substituting the contributions of the £'s in Eq. (llObl) . the contribution from 



the scattering interference within the QED theory can be written as the sum of two parts: 
a time-dependent scattering interference contribution 

aM«4/e i(£3d - e4/) '£i s i s £ 3 d4/{e- i(e3d - £4/) 5 - ^~ e *f^} (22a) 
+«I / a3de 4 ^- £M ) t £ lsls £ 4/ 34e- 1 ^-^^ - e l ^~ £ ^}, (22b) 
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which is non-zero, and a time-independent scattering interference contribution 

2£is ls [|a 3d | 2 £ 3(i 3 d + |tt4/| 2 £ 4 /4/], (23) 

which is identical to the one obtained within semiclassical theory. We rewrite the time- 
dependent interference contribution, Eq ( I22p . to the leading nonvanishing order in r as 

iCisis{BAf - e 3d )r [a&piifCaMfe***-**** - al^C^e^-^ 1 ] ± 0. (24) 

Therefore, for the particular combination of orbitals involved in the one-electron wavepacket 
and the orbital corresponding to the stationary electron in helium (see Fig. [T|), the time- 
dependent scattering interference between non-stationary and stationary orbitals is zero 
within the semiclassical theory and negligibly small in comparison to the total scattering 
signal in the QED theory. However, the time-independent scattering interference between the 
orbitals, i.e., the interference between Is and 3d orbitals and Is and 4/ orbitals, contributes 
equally to the total scattering signal in both the theories, as is evident from Eqs. ( 12 ip and 
( 123]) . The reason why the time- dependent scattering interference contribution to the total 
signal is so small in QED theory can be understood as follows: For probing the ultrafast 
motion, one has to satisfy Ae r C 1, where Ae = e^ — e^d is the characteristic energy 
scale of the electronic wavepacket, i.e., the pulse duration of the probe pulse should be 
much smaller than the characteristic timescale of the motion. Thus, the contribution from 
Eq. (|24|) . which is proportional to Ae r, is suppressed. 

On comparing the scattering contribution from the non-stationary electron with respect 
to total scattering signal, i.e., comparing Figs. |3^a) and EJ^a), one can easily distinguish 
the scattering contributions from the stationary electron and the interferences to the total 
scattering patterns in the QED theory. In Fig. [2]^a), the broadening of the scattering signal 
in the high Q region is the reflection of the contribution from the time- independent scattering 
interference between orbitals, whereas the wing type structures along the diagonal, which 
change as a function of the delay time, are a reflection of the contribution from the time- 
dependent scattering interference between orbitals, which is weak in comparison to the 
total scattering signal. On the other hand, on comparing the patterns shown in Figs. [21(b) 
and &b), one can only observe the broadening of scattering patterns, which reflects the 
contribution from the time-independent scattering interference between orbitals. 

Therefore, the patterns within the semiclassical theory are dominated by the scattering 
contribution from the time-independent interference between orbitals, and are not changing 
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FIG. 3: Scattering contributions from the non-stationary electron to the total scattering patterns 
in the Q x - Q z plane (Q y = 0) in helium. Scattering patterns obtained using (a) Eq. (|10c|) and (b) 
Eq. (jocj) at pump-probe delay times 0, T/4, T/2, 3T/4, and T, where the oscillation period of the 
electronic wavepacket is T = 6.25 fs. 

significantly as a function of the delay time. In contrast, the patterns within the QED 
theory are dominated by the scattering contributions from the non-stationary electron due 
to Compton scattering within the finite energy detection range of the detector. 

IV. CONCLUSIONS 



This work is devoted to understanding ultrafast time-resolved x-ray scattering from a 
sample containing a mixture of a non-stationary electron in the form of a one-electron 
wavepacket and one or more stationary electrons using the semiclassical theory and the QED 
theory of light-matter interaction. The contributions of the scattering interference between 
the non-stationary and the stationary electrons to the total time-dependent scattering signal 
are investigated in both the theories. Our investigations are based on the recent theory 
for time-resolved x-ray scattering to image the electronic wavepacket motion [40|. First, 
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we investigated different scattering contributions to the total scattering signal in both the 
theories and showed that the total signal can be decomposed into three main scattering 
contributions: first from the stationary electrons, second from the non-stationary electron 
and third from the interference between stationary and non-stationary electrons. In both the 
theories, the scattering contributions from the stationary electrons to the signal are identical, 
whereas scattering contributions from the non- stationary electron are completely different. 
In the QED theory, the scattering contributions from the interference depend on the energy 
resolution of the detector and the x-ray pulse duration. Therefore, in case of negligible 
energy resolution or extremely short pulses, QED theory provides identical contributions 
for the scattering interference as one obtaines in the semiclassical theory. On the other 
hand, if the pulse duration is not very short in comparison to the dynamical timescale of 
the motion and if the energy resolution is sufficiently high, the scattering interference in 
the QED theory does not provide identical result to the one obtaines in the semiclassical 
theory. It is important to note that most organic molecules and proteins contain mainly 
hydrogen, carbon, nitrogen and oxygen atoms and in such cases the scattering signal is 
dominated by the scattering from carbon, nitrogen and oxygen atoms. In these atoms, only 
two electrons are deeply bound core electrons, whereas other electrons are loosely bound 
valence electrons. When a pump pulse with broad bandwidth initiates excitation in such 
atoms it might be possible that more than one electron participates in the formation of the 
electronic wavepacket. In such situations, the time-dependent scattering signal would be 
dominated by the scattering contributions from the non-stationary electrons rather than the 
stationary and the interference contributions. 

Both the theories for light-matter interaction are illustrated by means of calculating the 
time-dependent scattering patterns for a one-electron wavepacket in helium. In helium, the 
pump pulse excites one of the electrons from the ground state and prepares an electronic 
wavepacket as a coherent superposition of the ls3d and ls4f eigenstates. The scattering pat- 
terns are computed for the non-stationary electron in the presence of a stationary electron. 
The time-dependent interference between the stationary and non-stationary electrons within 
the semiclassical theory is zero, and it is quite small in comparison to the total scattering sig- 
nal in the QED theory. However, the time-independent interference between the stationary 
and non-stationary electrons contributes identically to the total signal in both the theories. 
The patterns are dominated by the scattering contribution from the time-independent inter- 
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ference within the semiclassical theory, whereas the patterns are dominated by the scattering 
contributions from the non-stationary electron due to Compton scattering within the QED 
theory. Henceforth, the dynamical features of the patterns cannot be captured within semi- 
classical theory. We expect that our present analysis of TRI using ultrafast x-ray scattering 
will find several important applications for exploring ultrafast dynamics in nature. With 



the recent advent of novel light sources, we also believe that our findings will shed 
ultrafast electronic motion, for example, in atoms, molecules and biological systems 
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